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Abstract 



The ferromagnetic q-state Potts model on a square lattice is analyzed, for 
q > 4, through an elaborate version of the operatorial variational method. 
In the variational approach proposed in the paper, the duality relations are 
exactly satisfied, involving at a more fundamental level, a duality relationship 
between variational parameters. Besides some exact predictions, the approach 
is very effective in the numerical estimates over the whole range of temperature 
and can be systematically improved. 
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I. INTRODUCTION 



We consider the ferromagnetic g-state Potts model U on a square lattice, limiting our- 
selves to the isotropic, nearest-neighbour (n.n.) case, with no external field. The hamiltonian 
of the model is given by 

<ij> 

where Si — 1,2, ... ,q,i denotes the lattice sites, 6 SiSj is the Kronecker symbol and K = J / kT, 
with J > 0. We suppose that the lattice is made of m rows and n columns, and that free 
boundary conditions are applied. Our interest will be in the limits m, n — > oo. 

For q > 4 the model has a temperature driven first-order phase transition at K t = 
K t (q) = ln(l + y/q). At K t the free energy /, the internal energies u Q and Ud of the ordered 
and disordered phase, and the latent heat are known exactly Still at K t , more recently 
an analytic formula for the correlation length has been also derived |||4|]. 

Due to these exact results, the 2D Potts model has an important role in the study of 
the first order phase transitions 0. As a matter of fact its properties have been studied 
through several methods, such as analysis of high or low-temperature series-expansions || , 
resummation of large q expansions Monte Carlo simulations ||10-H . As a rule, the 



discontinuous character of the first order phase transitions makes the numerical calculations 
in these cases more uncertain, with respect to analogous calculations in the case of second 
order phase transitions. Moreover we have that some aspects of the discontinuous transitions 
are not completely understood from a general, theoretical point of view. 

On the other hand an accurate description of the Potts model transition, when the 
number of components q is very large, is provided by the Mean Field Theory [[Hifl. Indeed , 
it has been proved that Mean Field Theory is exact in the limit q — > +oo [[UJ. So it suggests 



itself that a reliable quantitative analysis of the properties of our model, at least for q > 4, 
could be also obtained directly from a proper extension of this theory. 

As a matter of fact, in this paper it is shown that considerable and systematic improve- 
ments of the Mean Field Theory can be obtained in the framework of the transfer matrix 
approach, through an elaborate version of the standard variational method for hermitian op- 
erators. This method is alternative to the variational Gibbs principle [[[9,0], which provides 
the analytical basis of the mean field procedure or, more generally, of the Cluster variational 
method [|TSj . 



The formal aspects of our approach are described in the next four sections. We will limit 
ourselves to q > 4. In the last section we present the results of some numerical calculations 
and we add some concluding remarks. 

II. OPERATORIAL VARIATIONAL METHOD 

Let Z mtn (K) be the partition function of our model. We are interested in the free energy 
per site, f(T), given by 

f(T) = -kT<j>(K), 
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where 

<f>(K)= lim lim — log Z mjn (K) (2) 

Now, if a = (si, s 2 , ■ ■ ■ , s m ) describes the spin configuration on a column, we can write Z m ^ n 
in the form 

Zm,n = XMH^U (*2)Lfo, ^3) •• ■ ^vOxM = (X, L^x), (3) 

<Tl,CT2,...,cr n 

where 

^-(cr) = e 2 ^<=i ^^+1 (4) 

and 

L(a,a')= X (<7)e KET=lSs »<x(<7') (5) 

is the symmetrized transfer matrix of our model. 

The matrix L(a, a') is hermitian. Therefore we can apply the spectral theorem and we 
can take into account of the extremal properties of its eigenvalues. We note also that, for 
our L(a,a'), the Perron- Frobenius (P.F.) theorem is valid [B5|. Furthermore we have that 
L(a,a') is positive definite, that is, for any vector £(<r), we obtain 

(£,LO = ^>)L(a,aW)>0. (6) 

ct,<t' 

In fact we can write 

m 

& m = /» ri( i + ^ - v^m^ 

a, a' i=l 

where 

f{a)=i{a) X {a). 

By expanding the product, we obtain a sum of 2 m terms, each giving a non negative contri- 
bution, due to e K — 1 > 0. 



Let \\{K; m) be the highest eigenvalue of L(a, a'). From (H) and (|3|) it follows that [19 



4>(K)= lim — log\i(K,m). (7) 
Now, by considering the functional given by the Rayleigh-Ritz (R.R.) quotient 
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we have 



sapF[£] = \ 1 (K]m), 
C 

the sup being attained for £(cr) = i>i{cr), where i[>i(cr) is the eigenvector of L(a,a') corre- 
sponding to Ai (K;m). So, for any vector £(cr), it results that 

lim - log %^ < <j>{K) (9) 

This inequality is the starting point of the operatorial variational method in Statistical Me- 
chanics pO| , pl| . In this approach, the principal (or ground state) eigenvector ipi (a) is usually 
described through an ansatz ip(a; a,/3, . . .), containing one or several unknown parameters 
which do not depend on m. Then, within the formal structure of the considered ansatz, the 
best estimate of f(T) is obtained by fixing the parameters through 

i- 1 i $,Lj>) 
sup hm — log — ~ ~ (10) 

a ,p,... m ^°° m (4>,4>) 

As a matter of fact, this approach, which is at the origin of the concept of the transfer 
matrix, is not really adequate to analyse the properties of a critical point. However the 
situation can be different in the case of first order transitions. 

In the following, fixing the attention on the 2D square Potts model, we will discuss a 
further elaboration of the above variational method. 

First of all, it is useful to bear in mind the meaning and the role of ipi{cr). In 
our lattice made of m rows, let us consider a strip S, of % adjacent columns, having n r 
columns at its right and n\ columns at its left. The reduced probability distribution (r.p.d.) 
Pi = Pj(o"i, 02, • • • , o"j) of a spin configuration on Sj , obtained by averaging over all the 
configurations at the right and the left of Sj, is given, in the limit n r ,ni — > +oo, by 

Df x ^i{vi)L{a 1 ,a 2 )L{a 2 ,az)...L{a i -i,(j i )'ilJi{a i ) 

PM, a 2 , . . . , a n ) = ' ( } 

for i > 2, while in the case of one column, we have the r.p.d. -Pi(er) through 

So, an ansatz ip{a) for V'i( cr ) has the role of a parametric description of the above averaging 
process leading to the r.p.d. Pj. 

Now, very reliable and simple descriptions, having however a different formal content, 
can be obtained by considering separately the high and low temperature regions. This fact 
drives us to introduce for ip(a) two formal structures, generally different, which we call 
ipdiv) an d i ) ( cr ) an d to consider ip(a) as made of two charts {il)^{a),ijj {a))^ associated to 
the ordered and disordered phase, which should match, if possible, at the transition point 
K t . The choice of the appropriate chart is demanded to the variational method. This point 
of view is a sligth generalization of the traditional approach, where a fixed formal structure 
for ip(a) is considered, and only the parameters are varied. However it leads quite naturally 
to variational approximations to <f>(K), in which duality [0,E3] is exactly satisfied. 
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III. ONE-PARAMETER ANSATZ'S 



In this section we introduce and discuss some one-parameter ansatz's of which we will 
make use in the variational approach based on equation fllQ|) . 

Let us fix the attention first on the disordered phase. For large T, due to the smallness 
of the correlation length, we expect that 

p ^-^r- (13) 

Then, with regard to the disordered phase, we introduce the simplest one parameter ansatz 
for ipi{o) in the form 

$\(t) = e^^S 1 ^,^!, (A(K) > 0) (14) 

that is, for K belonging to the disordered region, the result of the averaging process at 
the left and the right of a column, is represented simply by a renormalized ferromagnetic 
coupling A(K) between n.n. spins, such that A(K) ~ K for small K. 

This simplest ansatz can be improved in several ways. Since the averaging process pro- 
duces new interactions among the spins of a column, we need, in principle, other parameters 
related to these new couplings. For example, a first correction to ^ (cr) can be considered 
by taking into account also second-neighbour interactions. However, it is possible to stay 
with a one parameter ansatz and, at the same time, to improve $4(0), by applying the 
transfer matrix L to ^ (cr). As a matter of fact, if L is applied repeatedly to ^ (cr), we 
reach at the end the principal eigenvector. Since L is positive definite, we can consider, 
more generally, the action of th. So, due to the simple structure of V>]j (c), we will fix our 
attention on the sequence of one-parameter ansatz's 

$V), $V) = £*#V). ••• ,$\v) = LH<t\*) (15) 

The reliability of one approximate description of ifj± (er) can be tested by considering the next 
approximation in the sequence. 

Then we are led to the R.R. quotients 



JD (A K) _ (ffUff) (2) (A K) _ W,r4 d , 



which can be also written in the form 



>>u K) - an 



We note that 

Z^(A,K) = (^,L-^) (18) 

is, for v > 1, the partition function of a Potts(z/; A) model on a strip Yj v of m rows and v 
columns, having the same couplings as the original model ([!]), except on the first and last 
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column, where there is the vertical n.n. coupling (A(K) + K)/2. This model provides an 
approximate description P^ of the r.p.d. P u (o~i, c 2 , ... , o v ). 

We also remark that to Potts(y; A) is associated a q u x q v transfer matrix £ U (A, K), con- 
necting two adjacent rows of £„. We will call A^ } (A, if) the highest eigenvalue of £ U (A, K). 

In each R.R. quotient we can fix the parameter A(K), by considering the variational 
estimate (|T0|). We obtain then a sequence of estimates <p^\K) of 4>(K), given by 



^ (K) = sup lim — log 



1 Z^iAK) 
•"—*m Z]f {A. K) 



(19) 



sup log "[ M = sup ^ (A, K) {u = 1, 2, . . . ) 
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where the sup is taken over values of A which give ferromagnetic couplings, that is A > 
for v = 1, A > -K for v = 2, 3, . . . . 

It is useful to point out that, for < A, K < +oo, the P.F. theorem is valid for the 
reduced transfer matrices of finite order £ U (A,K), so that, in this range of the parameters, 
the Potts(z/; A) model has a unique ground state, as it happens in the thermodynamic limit 
for the original model, in the disordered phase. 

Now, let us consider the ordered phase. For finite m, even if large, ipi (a) is a well defined, 
unique, vector; furthermore P\{o~) has the same symmetries of the hamiltonian (JJ). However, 
in the limit m — >• oo, we have several phenomena as asymptotic degeneracy || of the highest 
eigenvalue of L, related to the occurrence of q ordered phases, macroscopic instabilities with 
respect to boundary perturbations or the existence of spontaneous long range order. So, an 
ansatz i/j (a) will be appropriate if these aspects can be correctly predicted. In particular, 
if we introduce reduced transfer matrices of finite order, as we did in the disordered case, we 
have to elude the P.F. theorem. 

We show that the above requirements can be satisfied through an approximate description 
not of the result of the averaging process leading to ipi(c), but of the process itself. 

Let us consider a column S. We replace the n r columns, with n r — ► oo, at the right of 
E, with one column E" of Potts spins a", erf, . . . , cr^, having an effective horizontal coupling 
H(H > 0) with the spins of E. Furthermore, adjacent spins erf are constrained to be in the 
same state, this last condition having the role to describe the ordered phase. The averaging 
process at the right of E is then simulated by the summation over all configurations on E". 
The same procedure is applied at the left of E, by introducing a column E'. So, we obtain 
the distribution 

P (a) ^^^n^^nKVle^S 1 ^, (20) 

O 1 a" 

where T is a normalization constant and 



m—l 



n(a>) = e" E -^n%^ +1 - ( 21 ) 



i=i 



6 



We are then led to the one-parameter ansatz ipo (cr) in the ordered phase, given by 

^)(a) = ^II(>>)j eT s Sl , Sl+1 ( 2 2) 

We note that P {&) is a r.p.d. obtained from the distribution 

P (a',a,a") = ^U(a' ,a)U(a" ,a)e K ^ s ^+^ (23) 

involving the strip S'USU Now, to P (cr' ', a, cr") is associated a c/ 3 x c/ 3 transfer matrix 
£ = Si, s'l \s' i+1 , Si + i, s'- +1 ), having many matrix elements equal to zero, due to the factors 

ni^i bs' v s' i+l i Y\T=i ^".s" i" ^ n * s a ^ ows to elude the P.F. theorem. In fact, the highest 
eigenvalue of i is q times degenerate; we have that £ is the direct sum of the q blocks made 
by the same q x q matrix, 

£(n) = £(si, s i+1 ; n) = £([M, s { , s i+1 , /i) (/i = 1, 2, . . . , q) 

and of the q(q — 1) blocks 

£(/i, //) = £(si, Sj+i; /i, fj,') = £(fi, Si, /i'l/i, Si+i, /-0 (//,//' = 1, 2, g // ^ //) 

made also by q x c/ coincident matrices, which differ from the previous ones. The highest 
eigenvalue of £ is contained in the blocks £{[x), which are associated to the q ordered phases, 
while £(/i, //') are related to the occurrence of a surface tension between different phases. 

The problem of selecting a particular £(/z) is physically equivalent to apply an external 
field hfj, to the spin state \i , on the boundaries £' and to take the limit m — > oo and then 
to consider the limit — > 0. Once this selection of a particular /i has been made, we see 
from ( pip that the horizontal effective coupling H gives rise to an external field 2H, applied 
on S to the spin state /i, obtaining then the long range order relative to the /i phase. So, we 
see that our parameter H has a role analogous to the effective field of the Mean Field Theory. 
By modulating through H the rigid boundary conditions on E' and provides the 

simplest description of the phase having a temperature dependent order parameter. 

As we did for the disordered phase, we can improve ^o\<j) through the action of L^. 
So, in the ordered region, we have the sequence of one-parameter ansatz's for ip\(<y) 

^\a) = L^\a), ... (24) 

and the relative R.R. quotients 

^■ K) --§0w> (v> - 1] (25) 

It is useful to point out that ijji 1 \<j) can be written in the form 



RETLi 1 s s>.s> , , + H T.T=i 5 s , ,s i +T Eili 1 s si<H+i 

( 7 ( ' 



v}P{(t) lim < 1 '''V. — ' m 



(26) 

lim ^\a;R), 
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and that 



(^U"- 1 ^) = Z { : +2 \H, K) (u > 1) (27) 

is the partition function of a Potts(z/ + 2; +oo, H) model on a strip Tj v+2 of m rows and 
v + 2 columns, with the same coupling as the original model (HD, except for the first and 
last column of sites, where we have a vertical n.n. coupling of infinite strength, and for 
the first and last column of rows where we have an horizontal n.n. coupling of strength H. 
This model provides an approximate description P„ of the r.p.d. P u ((Ti, <J2, ■ ■ ■ , oj) in the 
ordered phase. 

To the Potts(z/ + 2; +oo, H) is associated a q u+2 x q u+2 transfer matrix £„ + 2(+oo, H, K) 
having properties analogous to that of the previous £ = £ 3 (+oo, H, K). If we call 
\ < i +2 \+oo, H, K) its highest eigenvalue, we obtain the sequence of estimates 4>o(K) of 
4>(K) in the form 

^\K) = sup log S + 2 +00,I !rt\ = SUp ^\H : K) {v = 1, 2, . . . ) (28) 

By taking a vector ^\a) of the sequence (|l|) and a vector tpo\a) of the sequence (f2~H), 
we obtain an ansatz 

>m = ^fw)^\°)\ (29) 

with two charts, providing the estimates of 4>{K), given by 4>^\k) and 4>o\k) respectively. 
These two functions are defined for any K. For each K, the choice between the two descrip- 
tions is made by a further application of the variational method, that is by considering 

m^(4\K),^\K)) (30) 

We will say that two vectors ^\<j) and ^^(a) are exactly compatible if it happens that 
<j)j\K) and 4>o\k) match exactly at the transition point K t . In the next section, we will 
show that such vectors exist. 



IV. DUALITY 

A crucial property of the 2D Potts model is the duality relation which connects the 
partition function in the high and low-temperature regions. In terms of the free energy, it 
states that 

^K) = \ogx 2 + ^K*), (31) 

where 

e K -l 

x - 
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and K*, the point dual of K, is defined by 

K* 1 

e — 1 



The equation ( |3~1~D is an exact analytic constraint, allowing the location of the transition point 
K t (given by x — 1) on the basis of the non-analyticity of <fi(K) at this point 0. However, 
besides this formal aspect, we have no trace, in the above relation, of the mechanism leading 
to a transition from the disordered to the ordered phase. 

Now, we will see that further insights and additional aspects of duality are provided by 
the variational method, by means of the ansatz's introduced in the previous section. 

First of all we remark that, due to the fact that the ordered and disordered regions are 
described by a unique vector ipi{ a ) f° r finite m, the two charts of a proper ansatz should 
share some common property, as an expression of their common origin. As a matter of fact 
there is a property that allows to associate univocally tpo\a) with L^ip^\a) = ip^ +1 \cr), 
leading us to consider, as proper ansatz's with two charts, the sequence 

^)(a) = (^ 2) (a),^)(a)), ... = V), (32) 

and then the sequence of estimates <p^\K) for <f>(K), given by 

4>®{K) = max(^ +1) (K), 4>f{K)) (i = 1, 2, . . . ) (33) 

This property, shared by (a) and ^ (a), is an algebraic structure, as it will be shown 
in the following. 

We first consider ^ (a - ), for v > 2. We have seen that this vector defines the Potts(z/; A) 
model on a strip E„. Let us call r = (ti, t 2 , . . . , t v ) a spin configuration on a row of E„. The 
transfer matrix £ U (A, K), which connects two adjacent rows, can be written in the form || 

4(A K) = V V W U = q^a 2 x u -%(A, K), (34) 

with 

V V = {I + xU®)(I + xU^) ...(/ + xU^-V) [y > 2) 

_\ q a\I+W^){I+W^) {u = 2) 

V ~ \q^-\I + ![/«)(/ + IC/^-^C/ + if/F) ...(/+ i^ 3) ) (i/ > 3), (35) 

where 

A+K 

e 2 — l 
a = — 

and lljp , lfj) 2 \ . . . , U^" ^ are g v x q u matrices defined by 



Ul 2t - 1] (T,T') = ^=_l[5 t ^ (i = l,2,...,i/) 



C/f) (r, r') = ^ tii4i+1 J] 5^ (j = 1, 2, . . . , v - 1) 
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They satisfy the relations 

(U®) a = y/qU® (i = l,2,...,2i/-l); 
U®UP = UlpU® {\i-j\>2); 
UVug+QU® = U® (< = l ) 2,... J 2i/-2); 
U®U^*U® = U® (i = l,2,...,2i/-l) 

These relations define the Temperley-Lieb algebra generated by £/,; , {7^ , . . . , [/,; 1 , 
and then all the eigenvalues of £ U (A, K), with exception of their degeneracies. 

Now we come to "ipi (<j)(u = 1,2,...), which defines the Potts(i/ + 2; +oo, H) model. It 
is useful to consider this model as the limit, for R — ► +oo, of Potts(f + 2; R, H), obtained 
through the substitution of ^o\a) with (c; i?) of the equation (|26|) . Then, following the 
previous procedure, we have that to Potts(^; R, H)(v > 3) is associated a q v x q u transfer 
matrix £ U (R, H, K), which can be written in the form 

e v (R,H,K)=V v W v (R), (37) 

with 



V, 



(/ + hUi 2) ) (i + hU^- 2) ) (/ + xUi A) )...(/ + x£/r~ 4) ) (i/ > 4) 

(/ + hUP)(I + W7«) („ = 3) 

^ = e - 2 VpV~ 2 (i + + -U^-^I + -U^) ...(/+ -f^) (i/ > 3), 



where 



e H - 1 e R 

h - 



Taking the limit R — > +oo, we obtain 

4(+oo, if, if) = K^(+oo) = g^^~ 2 4(+oo, F, if) (39) 

with 

W v (+oo) = q^x v ~\l + -C/f) . . . (i + -f/f^) (40) 

So, the matrices t/jP and Uu 2u ~^ being suppressed by the limiting procedure, we remain 
only with the 2{y — 1) — 1 q u x q v matrices lll 2 \ llj) 3 \ . . . , u!) 2u ~ 3 \ lll) 2u 2 \ However, if we 
define 

Tjf = U$+V , (< = l,2,...,2(i/-l)-l) 

r(0 



we see from fl5E|), that these 2(z/ — 1) — 1 q u x q u matrices TO satisfy the same algebraic 



relations which are satisfied by the g 1 ' 1 x q u 1 matrices U^l^i = 1,2,..., 2(z/ — 1) — 1). We 
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only have with the U^ 1 a representation of different dimensionality of the same algebraic 
structure. Together with the structure of £ v (+oa, H, K) and £ v ^i(A, K), this is the common 
property which allows to couple the vector ^^(cx)^ > 3) with the vector $4 
leading then to the sequence (|32|). 

As a matter of fact, it follows from fl3"4|), (35), (38), fl3Tf ) and ( pif ) that the eigenvalues of 
4(+oo, H, K) coincide with those of £ u -i(A D (H, K*), K*), where A D (H, K*) is defined by 

A D (H,K*)+K* 

6 5 - 1 = \ (41) 

In particular, we have 

\{ u \+oo, H,K) = ^--x 2 ^\ { r l \A D {H,K*),K*) (v > 3), (42) 
V ^ 

which is a duality relation for the partition function of Potts(z/; +oo, if) and Potts(^ — 1; A). 

As a consequence, coming back to the variational method, we obtain the duality relation 
for the R.R. quotients 

\^ 3) (+oo,H,K) = x2 \ ( ; +2) (A D (H,K*),K*) 
\ { ; +2 \+oo,H,K) X X^ +1 \A D {H,K*),K*) 

Since, as H goes from zero to +oo, Ad(H, K*) goes from +oo to —K*, we deduce from (S3), 
(H), (H) that 

<j>^(K)= log x 2 + ^ +1) (K*), 

or, equivalent ly, 

0j +1) (K) = logx 2 + 0W(ir*), (44) 

the relative sup being attained at dual points, according to flHf ). So V>o (c) and ^"^(cr) 
are exactly compatible. 



The equation ( f44|) is the variational version of (23). Indeed, from fl44|) it follows that 

^(K) - 0i m) (K) = $ +1) (K*) - 4>f{K*) (45) 

Therefore, if at a point K we have, for example, (f)j +1 \K) > (f>y(K), then at the dual point 
it results that 4>o\k*) > <p^ +l (K*). So we conclude that the estimates <ffr\K) given by 
(B3I), satisfy exactly the duality relation (]3"ID, that is 



4> {i) (K) = log x 2 + 4> {i) (K*) (i = 1,2,...)- (46) 

We see that the mechanism of non-analyticity at K t can be described as a crossing phe- 
nomenon, involving (j)o(K) and <fij (K). 
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V. TWO-PARAMETERS ANSATZ'S 



The above results can be extended to two-parameters ansatz, which are useful to consider 
both from the computational and variational point of view. 

To fix our ideas, let us suppose that we need to go beyond the estimates (f)f'(K) and 
(j)j(K). We can consider then <j)^\K). The relative calculations require strips of four and 
five columns, to be compared with the case of three and four columns needed for (jr d (K). 
However we can also improve (p^' (K) without changing the number of columns and the size 
of the involved transfer matrices, but by adding a further parameter. In fact, by making use 
of the procedure of an approximate description of the averaging process, we can introduce 
the two-parameters ansatz for ipi(a) 

<& (3 V) = (^^^^^^^Je^ 1 W (47) 

a' 

with A,B > 0. 



o 9— 

A K 



o o o o 



o o- 



0- 



— — 

k Y b 



o o o o 



9 9- 



9 6- 



9 9 9 9 



9 9 9 9 



FIG. 1. The 3 and 4 columns Potts models originated by the two-parameters ansatz (f47|). They 
correspond to the partition functions zjp and leading to our estimate (f>^ of the free energy 
in the disordered phase. 

Then, by considering (ip d , Lip^ )> we are led to the estimate of 4>(K) 



<pf\K) = sn V limllog Z ^ B ^ , (48) 



where Z^(A, B, K){i = 3,4) are the partition functions of the Potts models shown in Fig. 
([!]). Therefore, we deduce from ( |T§D and (f£5|) 

4>f{K)<^\K)< <fi(K) (49) 
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We can improve analogously 4>^\k), which is determined through "ip^\a) = Lz'iffl (a), by 
introducing the two-parameter ansatz 

$?\a) = LH'P{a) (50) 

with 

$W(<r) = (j2 u (°'^)) e ^ EZ ~ llSa " St+1 ( G > °)> ( 51 ) 

cr' 

_____ ~' ( t l\ ~ f (2) ~'('2') ~' (1') 

where H{a',a) is given by (BTJ) . Through (?/v , Lipo )/(w >w )j we obtain the estimate 

= sup hm - log g^CgigLg) , (52) 

where Zs" (G, H, K)(i = 4, 5) are the partition functions of the Potts models shown in Fig. 
©■ 
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FIG. 2. The 4 and 5 columns Potts models originated by the two-parameters ansatz (|50|). They 
correspond to the partition functions and Z^ leading to our estimate (f)^ 1 of the free energy 
in the ordered phase. Double lines connect spins constrained to be in the same state. 

By making use again of the generators of the Temperley-Lieb algebra we have, as in ([[]]) . 

<pf\K) = log x 2 + ^ 2 \K*), (53) 

so that also the estimate of <f)(K), 

<f>'W(K) = max(tf 3) (ir),^ 2 )(ir)), (54) 

satisfies the duality relation ([31]). As a matter of fact, the parameters H and G are the dual 
expression of A and B, respectively. 

In principle, further improved two-parameters ansatz's can be obtained from the action 
of(Ls) 1 on the vector ^'^(a) = (&), (a)), leading to the sequence of estimates 
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VI. SOME NUMERICAL RESULTS AND COMMENTS 



Now we discuss the effective aspects of the variational method presented in the previous 
sections, by considering the estimates of lower order. 

We have calculated explicitly <f)f{K), <j>f{K), 4>'} 3) (K), (j) { o\K), 4> { o\K) and ff\K), 
obtaining then cj)^ (K) , (K) and <j)'^(K) } according to fl33|) and (0). We have also 
considered 4>^\A,K), which does not have a dual partner in our scheme, but, in any case, 
allows to obtain the simplest approximation to <f>(K), for every K. As a matter of fact, 
<fij(A, K) has a finite absolute maximum point A(K) for K < K t} while for K > K t the 
absolute maximum is taken at A — > +00. These points are also relative maximum points 
in a neighbourhood of K t , at its right and at its left, respectively. K t is an estimate of the 
transition point K t ; at K t we have a crossing between <fij(A,K) and 4>^\+oo, K ). 

The simplest estimate of our variational method is then given by 

4>V>{K) =max(0i 1) (A(K),K),0i 1) (+oo,ir)). (55) 
It turns out that <jyP>(K) satisfies quite well, even if approximately, the duality relation (j31~l), 




4.5 1 + -^-1- 

2 2.2 2.4 2.6 2.8 2 2.2 2.4 2.6 2.8 

K K 
(a) (b) 

FIG. 3. Our simplest estimate (full lines) of the free energy (a) and of the internal energy (b) 
at q = 100 are compared with the results from the Mean Field Theory (dashed lines). Also shown 
are the exact values at the transition point (small circles). 

when q is large. As a matter of fact, by making a comparison with the known exact results, 
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we have that <j)(°>(K) provides a quite reliable description, for every K, of the properties of 
the 2D square Potts model for q > 50, the description being more and more accurate as 
q increases. In this respect, our simplest approximation behaves like Mean Field Theory. 
However, from a quantitative point of view, we get a net improvement, since it turns out 
that the Mean Field estimates, which are exact in the limit q — > +00, differ significantly 
from the known exact results already at q of order 10 3 . Indeed, in the large q expansions, 
we have power series in 1/ y/q with very large coefficients 0, which are responsible of such 
deviations. In Figures (§)a and fl3|)b we plot, for q = 100, <p Q >(K) and the relative internal 
energy u^°\K), given by 



and we make a comparison with the analogous Mean Field quantities 4>mf{K) and umf(K); 
the small circles are the known exact values. 

From a computational point of view, in order to avoid q v x q v matrices of very high order, 
especially if q is large, we have made use of the relationship between the reduced transfer 
matrices and the matrices Uu , . . . , Uv 2v ~ 1 ^ which satisfy the relations fl36|) (see ([34]) , (35), 
(j37p and (38)). Now the relations (|36|) are also satisfied by the A u x A u matrices Uir given 
by § 

Qf 2 , • • • , a2v\a'n a' 2 , . . . , a' 2u ) = 5 aia / . . . 5 ai _ ia '._ i h a . a . +1 h a '. a '. +i 5 CH+2a '. +2 . . . 8 a2v a! 2v 

(56) 

(i — 1, 2, . . . , 2v — 1), where = ±1, a[ = ±1, 

h ++ = h = 0, = e~2, h_ + = e2 



A 



with 2 cosh A = q 1 ^ 2 . This is the 6- vertex representation of the Temperley-Lieb algebra. 

The eigenvalues of our reduced transfer matrices have been determined by making use 
of the representation Q55D, which allows to translate a Potts strip Sj into an ice-type (or 
6- vertex) strip and, as a by-product, to consider also non integer values of q. 

As we had to expect, the results of our calculations show that K t is the only crossing 
point of ^ +1 \K) with <$(K)(i = 1,2) and also of <p^\K) with <pT\K). Furthermore 
^ +1 \ K ) > <Po\K), <pf\ K ) > <Po 2) (K) for K < K t , so that 

^HSr *"»™-\z*w k>k, <57) 

In order to see the degree of accuracy which we obtain, we compare in Table | our results 
at K t , for several q, with the known exact values <fi(K t ) and —Ud(K t )/J. Of course, due to 
the duality relation, we reproduce exactly, through (j)^(K), (p^(K) and (p ^(K), the mean 
value (u d (K t ) +u (K t ))/2. 

We see that, for every q > 4, <p {0) (K t ), <p {1) (K t ), <p {2) (K t ) and ^ {2) {K t ) provide a 
monotonous sequence of very accurate estimates of <p(K t ), the accuracy being better as 




15 



TABLE I. Our predictions for the free energy <j>, the internal energy in the disorder phase 
and the specific heat in the ordered phase C Q calculated at the transition point Kt- The exact 
values for the free internal energy and the internal energy are also reported. 
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1 nn 
1UU 




9 /Ln^snfi 

Z.^UooUO 




3 909^94 




O. i OoooO 






2 407461 


2 898352 


3 204045 


3 429543 


3 758951 


4 803915 




2.408681 


2.899057 


3.204432 


3.429775 


3.759055 


4.803922 


0' (2) (^t) 


2.409378 


2.899363 


3.204563 


3.429841 


3.759078 


4.803923 




2.410306 


2.899597 


3.204631 


3.429867 


3.759084 


4.803923 


<f>(K t ) 


2.409849 


2.899522 


3.204615 


3.429862 


3.759083 


4.803923 




1.189595 


0.868464 


0.701564 


0.597935 


0.474256 


0.236934 


-u ( p{K t )/J 


1.249544 


0.911736 


0.728555 


0.615741 


0.483472 


0.237961 


-uf{K t )/J 


1.282729 


0.933513 


0.739933 


0.622203 


0.486151 


0.238115 


-u d (2) (K t )/J 


1.317667 


0.949882 


0.746210 


0.625067 


0.487035 


0.238139 


-u d (2 \Kt)/J 


1.413784 


0.971819 


0.751328 


0.626773 


0.487391 


0.238142 


-u d (K t )/J 


1.420754 


0.968203 


0.750492 


0.626529 


0.487353 


0.238142 


C?\K t ) 


4.5049 


5.4063 


4.6024 


3.8641 


2.9059 


1.1762 


d 2) (K t ) 


6.2350 


7.2927 


5.6956 


4.5144 


3.1923 


1.1969 


C {2) {K t ) 


9.9421 


9.8526 


6.6928 


4.9700 


3.3350 


1.2013 


c'o {2) {K t ) 


30.6324 


15.7009 


7.9572 


5.3728 


3.4165 


1.2022 



q increases. At q = 5, the errors of these estimates are of 0.25%, 0.10%, 0.05%, 0.02% re- 
spectively, while at q = 10, the sequence of the errors is 0.12%, 0.04%, 0.02%, 0.005%. As a 
matter of fact, by analysing the stability of our results for K ^ K t , we obtain that these 
errors are still rapidly decreasing, as we move from K t . 

For the internal energy, our predictions are characterized by larger errors, provided q is 
low and K is very near to K t . But our last estimate, through the two-parameters variational 
ansatz, is quite accurate for every K and every q > 15. At q = 20, — u d (K t )/J gives the 
exact value with an error of 0.1%. On the other hand, we obtain the errors of 7% and 2% 
at q = 5 and q = 10, respectively. However, as we move from K t , our estimates of the 
internal energy are again really quite accurate, even for low g's. This is shown in Fig.([|)a 
and Fig.(g)b, where we plot the functions (K) / J, -w (2) (K) / J and -w' (2) (K)/J. 

Through <p^\K) [i = 0,1,2) and (j)^(K) we have also determined the estimates 
C^'(K) and C ^ 2 '{K) of the specific heat C(K), which is not exactly known, even at K t . 
The sequence C^\K), C^'(K), C^ 2 \K), C'^ 2 \K) is monotonous and rapidly converging for 
every K, provided q > 30. For these g's, C w(K) provides a quite reliable estimate of C(K) 
for every K. This is also true for C^°\K), C^(K) and C^(K) as soon as q increases. For 
lower values of q, we have to move from K t in order to get the same accuracy. In Fig.([5D 
we plot C^(K),C^(K) and C'^{K) at q = 10, while in Table | we give, for several q, 
C { o '(K t ),C v '(K t ) and C Q K \K t ), that is the previous estimates at K t , evaluated in the or- 
dered phase (the predictions for the disordered phase can be obtained from these through 
duality) . 
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1.2 1.3 1.4 1.5 1.6 1 1.1 1.2 1.3 

K K 
(a) (b) 

FIG. 4. The monotonous sequence of our estimates —u^(K)/J, —u^ 2 \K)/J and —u^ 2 \K)/J 
of the internal energy at q = 10 (a) and q = 5 (b). The dashed line represents —u'^ 2 \K)/J, 
obtained from eq. (|5§|) . Also shown are the exact values at the transition point (small circles). 



On the other hand, we obtain a considerable improvement in our estimates of the internal 
energy and of the specific heat for lower q and K very near to K t , by accelerating the 
convergence of the already rapidly convergent sequence of our estimates of the free energy. 
Starting from <jy- x '(K), (p^(K) and ( 2 \K) and by making use of the Aitken A 2 method 
(or of the [1, 1] Pade approximant), we deduce the new extimate 

1 <t>( 2 )(K)-<j,W(K) 

We call u^\K) and C"^ 2 \K) the internal energy and the specific heat calculated through 
<f)"( 2 \K). The improvement which we obtain can be checked by comparing uJ 2 \K t ) with 

_ ///0 s ! ff (2) "(2) 

the known exact value. We give in Table I f {2 >(K t ), u d w (K t ) and C K \K t ) for several 



q. We see that — u d (K t )/J reproduces —Ud(K t )/J within the errors of 0.5% and Oa 
at q = 5 and q = 10, respectively, a net improvement with respect — w d (K t )/J. On the 
other hand, when u^\k) and C'^- 2 \K) are already quite reliable, <p"^(K) does not lead 

if (2\ 

to significant corrections. For q = 10 and q = 5, the function — u d (K)/J is plotted in 
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q=10 



20 

c 




FIG. 5. The monotonous sequence of our estimates C^\K), C^ 2 \K) and C '( 2 \K) of the spe- 
cific heat at q = 10. The dashed line represents C"( 2 \K), obtained from eg. (|58|) . 

Fig.([|)a and Fig.(f|)b, where its role is made quite evident. A stronger correction is obtained 
for specific heat, as shown in Fig.([|), where we plot C w(K), for q = 10. Our estimate of 
16 for C a (K t ) at q = 10, agrees essentially with the estimate of 18 obtained from large-g 
expansions |J, which is contradicted by the value 32 deduced from low temperature series 
||. For q = 15, 20, 30 the results of ref. || are confirmed quite well by our C (K t ). 

So, we see that the variational approach presented in this paper, besides providing exact 
predictions and physical or mathematical mechanisms for the first order transition of our 
model (for q > 4), is also very effective in the numerical predictions of unknown quantities. 
Our numerical analysis will be completed in a future paper and will concern properties as 
spontaneous magnetization and susceptibility, which are accessible through low temperature 
series expansions ||, but not through large-g expansions. On the other hand, since these 
two approaches give contradicting results in the case of specific heat, a further analysis of 
the estimates made in ref. would be appropriate. 

Furthermore we will explore more completely the consequences of duality, which in the 
variational approach proposed here appears at a more fundamental level. Through the 
reduced probability distributions and the variational parameters having dual partners, it is 
possible to study more completely the link between typical properties of the ordered and 
disordered phases. 

Our aim will be also to extend the above approach to the important case of the 3D Potts 
model. 
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